#'###########################################
#'                                          #
#'  ------------------------------------    #
#'  |     COIN     __                  |    #
#'  |            >(' )                 |    #
#'  |              )/                  |    #
#'  |             /  `----/            |    #
#'  |             \  ~=- /             |    #
#'  |   ~^~^~^~^~^~^~^~^~^~^~^~^~^~^   |    #
#'  ------------------------------------    #
#'                                          #
#'###########################################

rm(list=ls())

source("C:/UQAR/Recherche/Maitrise/Donnees/R/DataManagement/loadData.R")

STATIC_WORKING_DIR <- "C:/UQAR/Recherche/Maitrise/Donnees/Bylot/Fixe"

#'#
#'# Recuperation des donnees statiques et comparaison avec le gps
#'#


static.errors.raw <- loadData(STATIC_WORKING_DIR)


#'#################################################################################################
#'#                                                                                              ##
#'#                                 Proportion de localisations                                  ##
#'#                                                                                              ##
#'#################################################################################################

cols <- c("mistyrose2", "indianred3", "darkred")
ylim <- c(0, 85)

static.freqs.raw <- static.freqs(static.errors.raw)

freq.raw <- static.barchart(static.freqs.raw, col=cols, ylimits=ylim,
								main="R\u00e9partition des localisations par visibilit\u00e9 et classe de localisation")

print.multi(list(freq.raw))
#print.multi(list(freq.raw), save=TRUE, save.dest="C:/UQAR/Recherche/Maitrise/Donnees/Resultats/staticFreqs.png")

#save.png(freq.raw, "C:/UQAR/Recherche/Maitrise/Donnees/Resultats/staticFreqs.png")

#'#################################################################################################
#'#                                                                                              ##
#'#                                 Evaluation de l'erreur                                       ##
#'#                                                                                              ##
#'#################################################################################################


#'#
#'# Rayon d'erreur(68e centile) pour chaque visibilite
#'# 

colors <- c("mistyrose2", "indianred3", "darkred")

static.radius.raw <- static.radius(static.errors.raw)

radius.raw <- static.barchart(static.radius.raw, col=cols,
		main="Rayon d'erreur en fonction de la classe de localisation et de la visibilit\u00e9")

print.multi(list(radius.raw))
#print.multi(list(radius.raw), save=TRUE, save.dest="C:/UQAR/Recherche/Maitrise/Donnees/Resultats/staticRadius.png")

#'#
#'# Erreur globale statique en latitude et en longitude
#'#


static.error.long <- disp.boot(static.errors.filtered, stat.quantile,
		var=diffLong,
		quantileBy=quote(list(id, LC)),
		q=68,
		statBy=quote(list(LC)),
		key=c("LC"))
static.error.lat <- disp.boot(static.errors.filtered, stat.quantile,
		var=diffLat,
		quantileBy=quote(list(id, LC)),
		q=68,
		statBy=quote(list(LC)),
		key=c("LC"))

static.errLongLat <- rbindlist(list(static.error.long[, type := "Longitude"],
				static.error.lat[, type := "Latitude"]))

barchart(mean ~ LC, static.errLongLat, groups=as.factor(type), stack=FALSE,
		ly=static.errLongLat[,lb],
		uy=static.errLongLat[, ub], prepanel=prepanel.errorbar, panel=panel.errorbar,
		col=colors, box.ratio=3,
		xlab="Classes de localisation",
		ylab="Rayon d'erreur (\261 IC 95%)")

